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We study the head-on coUision of two equal mass, nonrotating black holes. Various initial 
' configurations are investigated, including holes which are initially surrounded by a common apparent 

, horizon to holes that are separated by about 20M, where M is the mass of a single black hole. We 

^\ • have extracted both £ = 2 and 1 = 4 gravitational waveforms resulting from the collision. The 

T-H ' normal modes of the final black hole dominate the spectrum in all cases studied. The total energy 

^jrj, radiated is computed using several independent methods, and is typically less than 0.002M. We also 

r- { ' discuss an analytic approach to estimate the total gravitational radiation emitted in the collision by 

^rfH I generalizing point particle dynamics to account for the finite size and internal dynamics of the two 

. black holes. The effects of the tidal deformations of the horizons are analysed using the membrane 

7—i ' paradigm of black holes. We find excellent agreement between the numerical results and the analytic 

estimates. 

y—^ • PACS numbers: 04.30. +x, 95.30.Sf, 04. 25. Dm 
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O ■ I. INTRODUCTION 

oo ■ 

o ■ 

I The spiralling coalescence of two black holes in orbit about one another is considered to be one of the most promising 
OS ' sources of gravitational waves 0| . The strong burst of gravitational waves resulting from such an event should be 

■ detectable by the next generation of gravitational wave detectors such as LIGO and VIRGO ||]. From observing these 
violent events in our universe, we expect to obtain important insights into astrophysics, gravitation and cosmology. 

^ ' In particular, such signals should provide the first direct and unambiguous evidence for the existence of black holes if 
the unique signature of the quasinormal modes |^ are excited. The information gained from the detected waveforms 

■ should allow one to reconstruct the astrophysical parameters of the system, such as the masses, spin, and orbital 
angular momentum and linear momentum of the colliding black holes, and the final black hole. Since LIGO and 
VIRGO are expected to begin taking data by the end of the decade, it is important to perform accurate calculations 

^ of the waveforms emitted during these events. Numerically generated waveform templates will be essential for the 
analysis of data collected by gravitational wave detectors. 

In a series of papers we investigate a special case of the black hole coalescence problem, namely the head-on 
collision of two black holes. On the one hand, the simplifying assumption of a head-on collision reduces the general 
three dimensional coalescence problem to a two dimensional axisymmetric problem, and is hence much more tractable. 
On the other hand, a head on collision can be regarded as an approximation to the last nonlinear stage of inspiralling 
coalescence-the final plunge. 

Our work extends and refines the earlier calculations of Dewitt, Cadez, Smarr, and Eppley ||7[-p^ (henceforth 
abbreviated as DCSE). Results from that collective body of work suggest that the normal modes of the final black 
hole resulting from the collision are excited and that the total energy released is typically less than 0.1% of the mass 
of the final black hole. However these numerical calculations proved to be very difhcult due to inherent coordinate 
singularities and numerical instabilities that plague the two black hole system. Also, the computer power available 
at that time did not permit highly resolved evolutions, and they did not have waveform extraction techniques at 
their disposal for determining gauge-invariant waveforms. For these reasons, DCSE quote their results as uncertain 
to within a factor of two, for example, in the total radiated energy [ p2[ . It is therefore imperative to revisit 
this important physical problem with the benefit of more powerful computers and improved analytic and numerical 
techniques developed over the intervening 15 years to calculate unambiguous waveforms and energy fluxes resulting 
from the collision. 
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Building on the work of DCSE, and more recent work involving distorted single black holes |14|-[16[, many of the 
numerical problems associated with colliding two black holes have been overcome in the present work. In particular, 
we have used a hybrid set of coordinates to resolve the axis and saddle point problems encountered by DCSE so that 
our evolutions are more accurate and more stable. The numerical code we have developed can evolve black holes 
with initial separation distances between ~ AM and ^ 20M, where M is the mass parameter defined to be half the 
ADM mass of the system. (We note that M is simply a convenient parameter to characterize the coordinates and 
only approximates the mass of a single black hole in the limit that the holes are widely separated.) We have applied 
more modern analysis techniques on the numerical data to extract waveforms and compute the total energies emitted. 
Using gauge-invariant waveform extraction we are now able to determine highly accurate £ = 2 waveforms and, for 
the first time £ = A waveforms. We also present analytic estimates of the total energy radiated in the collision by 
taking into account the finite size and internal dynamics of the two black holes. Comparisons of these results to 
the numerical solutions are valuable not just as a confirmation of our results, but more importantly, they provide a 
physical understanding of the numerical data. 

Sec. H provides an overview of the basic theoretical and computational methods used in this work. A more complete 
discussion of the initial data, the choices of coordinate systems, gauges, and the numerical methods we developed 
to handle the problems associated with the evolution of two black hole spacetimes can be found in a companion 
paper In Sec. Ill we present results from our numerical studies. First we investigate the highly nonlinear and 
dynamical near field region by looking at the evolutions of the metric data and apparent horizons. We further discuss 
the timing of the merging of the holes and its implications. Next we turn to the far field and discuss the extraction 
of gravitational waveforms and compute the total energy radiated from the collision process using various radiation 
indicators. In Sec. IV, we present a semi-analytic approach to estimate the total energy radiated and compare these 
results to the numerical calculations obtained in Sec. III. We conclude with Sec. ^ by summarizing our results and 
outlining a program to extend this work to more general black hole interactions. 



II. COMPUTATIONAL FRAMEWORK 

We use the 3+1 (ADM) formalism jl^ to write the Einstein equations as a first order (in time) set of differential 
equations for the dynamical variables "fij and Kij, the spatial 3-metric and extrinsic curvature respectively. The 
maximal slicing condition Tr i^T = is imposed throughout the evolution. Also, because the spacetimes we work with 
possess an axial Killing vector {d/dx^ = d/d(j)) all variables are independent of the azimuthal angle (j). 

This work, like that of DCSE, is based upon studying the axisymmetric evolution of the analytic Misner initial 
data [ p^ representing two equal mass black holes at the moment of time symmetry {Kij — 0). The spatial 3-metric 
for this data set can be written using cylindrical coordinates as 



df = [dp^ + dz^ + p^d(jP) , (1) 



where 



smh(n/i) \^ + + 2^j2 + _ 2^)2 j 

and Zn = coth(?T,/i). The free parameter n determines the total ADM mass of the spacetime and the proper distance 
between the two throats. The effect of increasing /i is to set the two black holes (centered along p = 0, at z = ± coth /i) 
further away from one another and to decrease the total mass of the system. 

The Misner data consists of two throats connecting two isometric sheets. The throats are spheres on which boundary 
conditions relating the metric across the two sheets may be imposed. Since the natural boundaries (the throats 
and a sphere surrounding the system far from the throats) do not lie along constant (z,p) coordinates, it is useful 
to introduce the body-fitted Cadez coordinates {rj^£,) where rj is the logarithmic "radial" coordinate and ^ the 
"angular" coordinate. The advantage afforded by this set of coordinates (shown in Fig. |^) is that they are spherical 
near the throats of the black holes and also far away in the wave zone, thus allowing us to deal with throat boundaries 
and asymptotic wave form extractions in a convenient way. The disadvantage is that the coordinate transformation 
introduces a singular saddle point at the origin [z = p = Q) that is not present in cylindrical coordinates. This creates 
certain numerical difficulties that require special treatment as we discuss briefly in the remainder of this section. 
The success of our methods depends critically on utilizing both sets of coordinate systems (cylindrical and Cadez) to 
advantage. 

We have investigated a number of different numerical schemes to solve the problem of colliding two black holes 
head-on. The basic idea that evolved from our investigations is to solve for the Cadez metric components, which we 
write as 
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7»j =* =^ { C B ] (3) 
D sin^ ^ / 

in the coordinate order (77, ^, 0), on the Cadez grid and use a shift vector to set C = dtC = 0. (We note that 
only the conformal metric components are evolved. The conformal factor 5''* = ^'j^^/ J, where J is the Jacobian of 
the Cadez/cylindrical coordinate transformation, remains constant in time.) This choice for the shift vector has the 
advantage of a diagonal 3-metric which helps to suppress the axis instability and simplifies the equations of evolution 
and the extraction of invariant gravitational waves in the far field. Furthermore, with this approach it is possible to 
define variables for the two black hole system that obey the same evolution equations with similar boundary conditions 
as the single distorted black hole code developed in previous work In fact, the two black hole code in its final 

incarnation evolved from the code we developed for distorted axisymmctric single black hole spacetimes and much of 
the discussion in is directly applicable here. 

The difficulty with Cadez coordinates is the singular saddle point located within the computational domain at 
the origin z = p = (see Fig. |l|). We evolve data near the saddle point by taking advantage of the fact that the 
spacetime metric components in cylindrical coordinates are smooth everywhere, including the saddle point. We can 
therefore define a cylindrical coordinate "patch" to evolve the cylindrical metric and extrinsic curvature components 
on the Cadez grid over regions near the saddle point. The two sets of components, Cadez and cylindrical, are evolved 
everywhere independently of each other (except for the coupling at the patch boundaries) on a single Cadez grid. The 
nonsingular cylindrical components are then used to correct the singular Cadez components in the patched region 
using the general tensor relations T/^- = {dx^ / dx'^){dx^ / dx'^)Tki. The Cadez components, in turn, provide corrections 
to their cylindrical counterparts everywhere else, helping to suppress the axis instability that is inherently present in 
the cylindrical coordinate system possessing a nondiagonal metric. A more detailed discussion of this procedure can 
be found in 

Our code was subjected to a number of tests, including matching waveforms to perturbation theory as we do in 



Sec. pi B l| . We also performed various convergence studies using 100 (27), 200 (35) and 300 (55) radial (angular) 
zones. We have shown that the convergence rate for the total radiated energy is quadratic in the grid spacing and, 
more specifically, differences in the dominant £ = 2 waveforms between the 200 and 300 radial zone evolutions is on 
the order of just a few percent. We discuss the accuracy and reliability of our calculations throughout this article 
when appropriate and refer the reader to Ref. p] for more details. 



III. NUMERICAL RESULTS 



We have investigated six separate cases of the Misner two black hole data sets corresponding to different values of /i. 
The physical attributes of the initial data for these six cases are summarized in Table | where we show M — Madm l'^, 
[ie., half the ADM mass), the proper distance between the two throats, and whether the data contains a single global 
apparent or event horizon surrounding the two holes. In this section, we categorically discuss the dynamics of both 
the near and far field regions for the various cases. 



A. Near Field 



1, Spatial Metric 



First we present results for the case /i — 1.2, which is a data set corresponding to two holes that have already 
merged initially. The initial data contains an apparent horizon that encircles both throats. In Fig. ^ we show the 
conformal metric function 7^^ = A ai the coordinate time t = 2bM . Notice that a sharp peak surrounding the hole 
is developing. This peak develops essentially spherically around the two throats from early in the evolution, showing 
that the system behaves as a single black hole from the outset. These results are similar to those observed in studies 
of single throat spacetimes p^ , p^ . The reasons behind the development of the peak is clear: As the coordinates are 
dragged into the hole, the proper distance between radial grid points increases rapidly towards the throat. However, 
as shown in Fig. ^ as the lapse goes to zero in the region near the throats, "freezing" all motions there. Hence we see 
only the growth of the proper distance between grids in coordinate time in the region near the horizon, developing 
a sharp peak. Such grid stretching effects present one of the main difficulties in evolving black hole spacetimes in 
any numerical simulation utilizing a singularity avoiding time slicing. During the course of evolution, the dramatic 
change in the radial metric function becomes increasingly more difficult to resolve numerically, and the region of 
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causal dependence eventually becomes much smaller than the coordinate grid spacing (as the radial metric function 
increases and the lapse decreases there). 

For comparison, we show in Fig. ^ the lapse for the case /i = 3.25 where the throats are much more separated, 
forming two black holes initially. The display corresponds to an early time {t — 22. 5A/) in the evolution of the system 
where the two holes are acting essentially independently of each other as they begin to fall together, and the collapse 
of the lapse locally around the throat reflects this fact. After the holes begin to coalesce, the lapse collapses spherically 
around both throats which are contained within the final black hole. 



2. Horizons 



We have already observed that at late times the intrinsic geometry of the apparent horizons oscillates at the normal 
mode frequency of the final black hole. This aspect of apparent horizon dynamics is discussed in detail in Ref. 1^ for 
three different spacetimcs including colliding black holes. Here we look at horizons in the context of distinguishing 
one from two black hole data sets and estimating the mass energy of radiated gravitational waves. 

One can define an effective mass of a black hole based on its apparent horizon via the relation 



where Ah is the intrinsic area of the apparent horizon. This relationship gives a lower limit for the mass of the black 
hole, since we know the apparent horizon should lie within or (in the case of stationary black holes) coincide with the 
actual event horizon. When the system begins to settle down Eq. (^ provides a good estimate of the black hole mass 
as the apparent horizon will lie very close to the event horizon. 

In Fig. II we show the evolution of the horizon mass computed from Eq. (^ for two different cases. The numerical 
data for the case /i = 1.2 is shown as a solid line. The mass of the hole Mh is normalized to units of the total 
ADM mass of the spacetime, so ideally we should have Mh < 1 for all time. However, because the horizon is always 
found near the peak of 7,,^, is extremely sensitive to the precise position of the located horizon as discussed, for 
example, in Ref. |^ . Small errors in the height and shape of the metric functions arising from inadequate resolution 
can exaggerate errors in the area of the horizon derived from these metric functions. A typical effect of this problem 
is to overestimate the horizon mass Af/,, after the metric functions become too sharp for the grid to resolve. (This 
effect has also been discussed in Ref. po| where an apparent horizon boundary condition was used to circumvent this 
problem.) 

Also shown in Fig. ^jis the result for the case /z = 2.2. In this case the initial horizons are the two distinct throats 
and contain only about 79% of the total mass of the spacetime. Because we use a slicing condition in which the 
lapse is zero on both throats, each throat remains a marginally trapped surface throughout the evolution. We track 
this surface only until another trapped surface forms across the equator (z = 0) to surround both throats. This new 
horizon, when it forms, contains essentially all the mass of the spacetime except for a small amount carried by the 
radiation. Ideally one should see the horizon mass leveling off just below the total ADM mass of the system, with 
the difference being accounted for by the energy carried away by gravitational radiation. However the errors in the 
horizon mass, due to effects discussed above, are large enough to hide this small amount of radiation which is less 
than 0.1% of the total mass. Fig. || also shows a spurious feature appearing at the time of the formation of the new 
horizon. This feature is due to numerical difficulties with locking on to the new horizon which first appears near the 
coordinate singularity present in the Cadez coordinate system. After a brief period, the system settles to the final 
black hole configuration (apart from the numerically induced slow growth in the horizon mass as discussed above.) 

The results for holes with wider initial separations are similar to the case fi = 2.2, except that the time scale for 
the merging of the holes is longer. Thus, by the time the holes have merged, as measured by the appearance of an 
outer horizon, the horizon mass is overestimated by larger amounts. In the case /i ~ 3.0 this effect is about 20%. 
Note that the wavelength of the quasinormal mode of a black hole scales as its mass, so this is consistent with the 



results of Sec. (HIB 1) which show that the extracted waveform for these cases have wavelengths slightly longer than 
expected for a final hole of mass Madm- 

Although apparent horizons can confirm when two black holes have definitely merged, the presence of two distinct 
apparent horizons does not guarantee that the holes are separate, as a global event horizon may still surround the 
individual holes. The claim that the larger /i cases represent two distinct black holes is supported by computing the 
area of a 2-sphere defined on a constant 77 surface that just encircles the saddle point and therefore both holes. In 
Fig. ^ we plot the area of this surface on the initial time slice as a function of separation (and mass) parameter /i. The 
horizontal line is IGtt, the area corresponding to a surface representing an effective mass equal to the ADM mass (our 
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units are such that the ADM mass is normahzed to unity). The point of intersection is /i ~ 2.1, and Fig. ^indicates 
that for larger values of ji the data represents two black holes. 

A more convincing and precise argument can be realized by integrating light rays radially outward along the equator 
starting from midway between the black holes. If the photons "escape" , then clearly the system contains two disjoint 
event horizons initially and hence two separate black holes. (Whether photons escape or not is determined by their 
position relative to the apparent horizon at late times, typically ~ 80Af . At such late times the apparent horizon 
is expected to lie near the event horizon.) We have performed several such studies in an attempt to find a critical 
value for that separates the initial data into one or two black hole sets. Our results indicate that this value is 
approximately 1.8. For /i > 1.8 the photons escape, while for < 1.8 they do not. Furthermore, we have recently 
developed a method for tracing out the actual event horizon surface [pT| and found results consistent with those 
obtained by integrating photons. 



In Fig. [7| we compare different timings of the collision for each evolution considered. The meaning of the various 
timings is as follows: tnewt is the Newtonian free fall time for two point particles to collide from rest and separated 
initially by a distance L/M ^ thor is the time at which the first merged apparent horizon appears, and trmg is the time 
at which the first (negative) peak in the £ = 2 gauge-invariant waveform reaches the detector at r = 40M, indicating 
the (retarded) onset of the quasinormal ringing of the final black hole. First we note the remarkable coincidence 
between the free-fall times to collision and the time required for the apparent horizons to merge. This agreement 
between such different indicators of the coalescence time is very satisfying. 

Another interesting feature is the timing of the onset of the quasinormal ringing. Note that for ^ below 2.2 
{L/M < 8.9) the ringing begins at about the same time regardless of the separation between the holes. For holes 
separated by more than this, a delay in the onset of ringing becomes apparent. The fact that the ringing begins at 
the same time regardless of their separation is not surprising for below 1.8. As we discussed above, the two throats 
are really just a single distorted black hole and the separation between throats is physically irrelevant. For fj, between 
1.8 and 2.2, as the two holes are initially disjoint, one might expect a delay in the onset of ringing. However, the 
important physical property of the system governing the quasinormal ringing is the gravitational scattering potential 
barrier surrounding the holes, not the position of the event horizon itself. Since the peak of this potential barrier 
is located near r — 3M for a Schwarzschild hole and the horizon is at r = 2M, we expect the potential barrier of 
each hole to merge into a single one before the horizons do. Therefore, for black holes that are initially close enough 
that their potential barriers have effectively merged, we expect to see a system that behaves essentially like a single 
distorted black hole in terms of quasinormal mode ringing, even if the event horizons are distinct. From Fig. ^ it is 
clear that this transition from merged potentials to distinct potentials takes place near fi 2.2. 



The main method we use to calculate waveforms is based on the gauge invariant extraction technique developed by 
Abrahams and Evans [|l^ and applied in Ref. ||l^ to black hole spacetimes. The basic idea is to split the spacetime 
metric into a spherically symmetric (static) background and a small perturbation in the region where the curvature is 
dominated by the mass content of a small compact object. We first expand the metric perturbation in m = spherical 
harmonics Yio{9) and their tensor generalizations. The Regge- Wheeler perturbation functions are then extracted from 
the numerically computed metric components and used to construct the gauge invariant Zerilli function ip. (See ]T3t for 
a detailed discussion of this procedure.) ip represents the wavelike part of the metric that is radiative at large distances 
from the source and is commonly used in semi-analytic calculations of black hole normal mode frequencies Q . The 
asymptotic energy flux carried by gravitational waves can be computed from 



independently for each i mode contribution for the normalization we use for ip. 

For all of the cases studied in this paper we have extracted both the £ — 2 and £ = 4 waveforms at radii of 
30,40,50,60, and TOM. (Coordinate positions corresponding to physical distances in units of M are approximated 



3. Collision Timings 



B. FAR FIELD 



1. Gauge-Invariant Waveform Extraction 




(5) 
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from the initial data in the asymptotically spherical far field as r = By comparing results at each 

of these radii we are able to check the propagation of waves and the consistency of our energy calculations. 

In Fig. H we show the £ = 2 waveform (solid line) extracted at a radius of 40M for the case /i = 1.2. This result 
is similar to waveforms extracted from simulations of single distorted black holes (see, for example, Ref. jl^) even 
though the initial data sets and coordinate systems used are significantly different. Of course since there is a horizon 
surrounding both throats, we expect this system to evolve as a single black hole from the outside. As a perturbed 
single black hole system, we also expect the quasinormal modes of the black hole to be excited. The dotted line in 
Fig. ^ shows the fit of the lowest two (fundamental and first overtone) £ = 2 modes of a black hole of mass (2M), over 
the range 70 < t/M < 160, obtained from Refs. |^,^. The fit is excellent, showing that the normal mode is the 
dominant part of the emitted radiation. We note that the first overtone quasi-normal mode is much more strongly 
damped than the fundamental, and hence does not contribute appreciably to the fit at late times. Its main effect is 
to increase the accuracy of the fit to the first peak in the extracted waveform. 

Next we discuss the case /j, = 2.2, for which there are no initial common apparent nor event horizons. In Fig. ^we 
show the £ — 2 extracted waveform for this case. The solid line shows the waveform detected at a distance r = 40M 
and the long dashed line shows the waveform extracted at r = 60M. The wave is clearly propagating away from the 
hole at light speed with essentially invariant shape and amplitude, with a wavelength of 2 x 16. 8A/, confirming the 
original findings of Smarr and Eppley [l^ ]. However, our more accurate code now allows us to go beyond estimating 
the wavelength and to fit quantitatively the waveform to results known from black hole perturbation theory. The 
short dashed line shows the result of fitting the r — AOM waveform (in the range GAM < t < 160Af) to a linear 
combination of the fundamental and first overtone of the 1 = 2 quasinormal mode for the final black hole with mass 
2M. The fit is quite good, matching both the wavelength and damping time, showing that the final black hole mass 
is indeed very close to the total mass of the spacetime. 

In Fig. ^ we show the more difficult i = A waveform for the same case = 2.2, extracted at the same radius 
r = AOM. Again this waveform has been fit (over a similar range) to a superposition of the fundamental and first 
overtone £ = A quasinormal modes of the black hole. Although the fit to the extracted waveform is reasonably good 
in terms of the wavelength and damping time, this waveform is rather sensitive to the computational parameters such 
as grid resolution and the extent of the numerical "patch" of cylindrical metric functions covering the saddle point. 
As discussed in Ref. Q, the amplitude of the £ = A waveform can vary by about a factor two over a wide range of 
patch, diffusion, and resolution parameters. Future refinements of this code may allow us to make more definitive 
predictions of this difficult waveform extraction. 

Finally, in Figs. [l| and |l| we show the £ = 2 and £ = A waveforms respectively for the case = 2.7, where the holes 
are well separated by about 12. 6M initially. The solid lines are the waveforms extracted at a distance r = AOM. In 
this case the fits to perturbation theory are still reasonably good, but are not as close as the calculations performed 
for holes that are initially closer together. The wavelengths of the extracted waveforms are somewhat too long, and 
this can be understood as a numerical artifact of our methods. The calculation must run for a longer period of time 
before the onset of quasinormal ringing, so the peak in the radial metric function becomes more difficult to resolve. 
This leads to an error in the longitudinal (spherical) part of the field, that causes the effective gravitational scattering 
potential to be somewhat different from the true potential. (See also Ref. |16 for a discussion of this point.) Since 
this potential is critical in determining the quasinormal frequencies of the system, the normal modes are generated at 
slightly different frequencies. 



2. Other Radiation Indicators 



The Newman-Penrose scalar 



*4 = Raf}^sl°'m^k''rrf 



(6) 



provides another approach that can be used to treat the problem of radiation extraction. The vectors k and I are 
orthogonal real vectors defined by adding and subtracting a spacelike unit vector with a unit timelike vector. The 
vector TO and its complex conjugate to are orthogonal null vectors tangent to the surface of a 2-sphere representing 
the wavefront of an outgoing shell of radiation. The basis set we have chosen to work with is the following 
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-a + '^'^^/lp'^, ■^'^^/A, 0, 
-a--^^^/A(3'^, -^^^/A, 0, 



(7) 
(8) 
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(9) 
(10) 



Far from the source ^'4 represents an outward propagating wave and is therefore naturaUy normal to a 2-sphere of 
constant "radius" rj. The total radiated energy loss can therefore be estimated by pa 



dE _ 1 



(11) 



where the integration is over a 2-sphere (with an area of 47rr^ and a surface element df2, where r ^ ^J^f^/M 
*2/B/Af) lying in a spacelike hypersurface surrounding the radiating system. 

A third method that we have used to track gravitational radiation is based on the Bel-Robinson vector 



(12) 



where Eap and Bap are the "electric" and "magnetic" components of the four dimensional Riemann tensor. 

is the 3-dimensional Levi-Civita permutation tensor and n" is the unit vector normal to the space- 



like t = constant hypersurfaces. Although p'~' is constructed in a manner formally similar to the Poynting vector of 
electromagnetism, we note that it is not a physical momentum vector for gravitational waves, since its units differ 
by M^. Nevertheless, has been proven to be effective in q ualitat ively tracking gravitational radiation |11 27tl, and 
as demonstrated by DCSE in Ref. | 
good quantitative agreement with t 



12[ | and again in section [II B 3 , the radiated energies computed using 
rose computed using other radiation indicators. Since 



are m 

is dimensionally a flux 



vector quantity, the energy loss can be approximated by integrating the radial energy flux over a closed 2-sphere 
as described above for ^'4 



dE 



1 

47r 



dt' 



'(±y'b^r,|/2 



H 2 



r^dn. 



(13) 



The choice of sign in the integral of (|T^) is taken to be the sign of \E'4. Equation (|T3|) is motivated by the asymptotic 
form p^r-y 2|vl'4p applicable to the case of monochromatic waves in linearized Minkowski spacetimes This 
construction yields results that are consistent with the integrals of and ^'4 in the asymptotic far field. 

It is informative to compare the waveforms obtained by the three different methods. In Fig. ^ we show the 1 = 2 
Zerilli function (solid line) , ^'4 (dashed line) and p"^ (dotted line) at a fixed point r = TOM along the equator for the 
case /i = 2.2. ^'4 and p^ have been normalized to the same scale of ip by matching the amplitudes of their maximum 
peaks. The £ = 2 fundamental quasinormal mode is clearly present and dominant in all three signals. 



3. Energies Radiated 

The total radiated energy E can be computed from the Zerilli function using Eq. (|^). We display these results 
in Fig. |lj. The six clusters of unconnected symbols represent the six numerical simulations corresponding to the 
different parameter values. Each of the five symbols within a cluster corresponds uniquely to the total integrated 
(. = 2 energy computed at the five different wave detectors. For reference, the early results of Smarr and Eppley are 
plotted as large crosses with error bars suggested by Smarr |Q. Within the large errors quoted, those early results 
are remarkably consistent with our more accurate results. 

Clearly the results in Fig. ^ show two distinct regimes, as denoted by th e arrow s in the upper part of the figure. 



For /i < 1.8 the initial data contains one black hole, as discussed in section [II A 2, and the energy radiated falls off 
exponentially. We have fit the energy output to an exponential for /i < 1.8 and find that the results are approximated 
well by the formula E = 3.13 x 10~^ exp(4.852/i). For /i > 1.8 there are two holes and the energy radiated is 
somewhat independent of the initial separation. Three lines based on analytic and semi-analytic calculations that 
treat the system as two black holes are also shown, as we discuss in more detail in the next section. The curve labelled 
"DRPP calculation" is based on Ref. , the result labelled "Reduced Mass Correction" takes into account the finite 
mass of a black hole falling into its partner, and the dashed line is the complete semi-analytic calculation discussed 



in Sec. IV, accounting for other effects. Finally we note that the energy radiated is very small compared with the 
upper limits based on the horizon area theorem (represented by connected circles in Fig. |lj), as discussed in the next 
section. 
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Table |l^ compares the total radiated energy computed using three different radiation variables: the Zerilli function ip, 
the Newman-Penrose scalar \l/4 and the Bel- Robinson vector . The energies in Table || are results from simulations 
resolved with a 200x35 grid and are normalized to the total ADM mass of the spacetime (2M). The range of values 
quoted for the numerical calculations for each of the six cases represent the range of data across the different detector 
locations, with the lower (upper) limits corresponding to the outer (inner) detectors. As a whole, the results are 
remarkably consistent. The various numerical constructions differ significantly only at the innermost detectors for 
the low fi cases. These deviations are attributed to several effects. First, the near zones are characterized by stronger 
highly distorted behavior than the asymptotic wave zones, making this region more difficult to resolve accurately. 
and are curvature type variables that have an explicit dependence on first and second order gradients of the metric 
components. Hence, they are more susceptible to numerical inaccuracies than the Zerilli function which depends 
only on first derivatives of the 3-metric. Secondly, the curvature quantities are projected onto a coordinate based 
tetrad, and the Cadez coordinates are distorted from sphericity in the inner regions of the grid, obscuring the physical 
interpretation of these quantities. 



IV. ANALYTIC ESTIMATE 



To gain a physical understanding (and confirmation) of the numerical results, we outline a procedure to estimate 
analytically the total radiated energy. Our approach is based on the well-studied problem of a test point particle 
originally at rest at infinity plunging into a Schwarzschild black hole ^9 28 30|-|3^. For the test point particle prob- 



lem, [ |29| combined a Newtonian quadrupole moment calculation with the linearized theory of Landau-Lifshitz |33| ] to 
find the in-flight radiation 

^ 105 \ m) ' ^ ' 

for infall from cx3 to r = 2M, where m is the mass of the test point particle and M the mass of the black hole with 
m <ti M. In 1^ the total radiation for the same test point particle problem is obtained using black hole perturbation 
theory (see e.g., [|| and references therein). The result is comparable psf 

? 

£; ==0.0104—. (15) 
M 

Our approach is to adopt the general relativistic result ( [Tsl ) and modify it to include correction factors so that it can 
describe the two black hole collision. In the following, we shall discuss correction factors due to (A) m is not much 
smaller than M, (B) the infall is not from infinity, and (C) the black hole, unlike a point particle, has a finite size 
and internal dynamics. 



A. mass scaling 

Before we go into the various correction factors, it is useful to understand why E in (^5|) is proportional to m? /M . 
For m M , the quadrupole moment of the system is / ~ mr^ where r is the radial distance between m and M . 
The gravitational wave luminosity is given by /, the third time derivative of /. In the Newtonian approximation, 
f - y/2M/^ and r - M/r^, we have 

Lcc'P r^m^{rrf r^m^(J^^ . (16) 

The total energy radiated is 

E = / LdtK -^strong field X iJ^strong field • (17) 



The integral is evaluated at the strong field region, as most energy is released towards the end when m is falling near 
the horizon of M. Putting r = 2M into (|l6|) for L in the strong field region, and (5tstrong field ~ M, we have 

9 

TJl 

Eoc-, (18) 
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as in (|5|). 

When m is not much smaller than Ai , it is more accurate to use 

/ - Air^ (19) 

where ^ = (mM) / (M + m) is the reduced mass of the system. Hence the considerations above (equations ( [l6| ) to ([l^)) 
suggest that Eq. (15) should be changed to 



E = 0.0104- 



M 



(20) 



as the gravitational wave energy output for the case when m is not necessarily much smaller than AI. Notice that for 
m = M, the /i^ of Eq. (|20|) introduces a quite significant factor of 1/4. 



B. finite infall 



The simple quadrupole approximation of (|l^) also suggests how the expression should be modified when the infall 
is not initially from infinity. Starting from rest at a finite distance Tq reduces the velocity f, which enters (^6|) as r^. 
Denote the energy radiated in such a case as E^^ , 



Eq. i ^dj ) is hence modified to be 



f.2M 

Er, ^ / Ldt^ / {L/f)dr . (21) 



E = Fr X 0.0104^ (22) 

Er Sr^^ '''{^Yd^ 



E„ 



loo '^{''^ydr 



with 



2M 1 /2M 2M 



1-2M 



(24) 



Since can be as small as a few M in the numerical simulation, we must use the relativistic expression for r 
in Schwarzschild coordinates in Eq. (p^). Equations (|3|) and (|4|) represent one way of extending the quadrupole 
formula to the highly relativistic regime. There is no unique way to do the extension, we have just picked a way 
convenient for our present purpose. This correction factor Fr^ represents two effects: (i) there is less time to radiate 
when falling from a finite distance, and (ii) the infalling velocity is smaller. The latter effect is much more important. 
In Fig. |l^ we plot Fr^ vs. Tq, the initial Schwarzschild coordinate of m, covering the range of ro used in the numerical 
simulation. 



C. internal dynamics 

All considerations up to this point are the same, independent of whether the infalling object is a point particle 
or a black hole. In the following, we consider correction factors due to this difference. Although at the end we will 
extrapolate to two black holes of equal mass m — AI , for both generality and convenience of discussion, we think of 
the situation as a hole with mass m falling into a hole with M > m. 

As far as the gravitational wave output is concerned, the most important difference between a point mass and a 
black hole is that a black hole has internal dynamics. There are more channels that the initial gravitational potential 
energy in the system can dissipate into. Such dissipations decrease the kinetic energy and hence the velocity of the 
infalling hole. Hence fewer gravitational waves are generated. There are various mechanisms causing dissipation, 
which we shall describe separately. Of course, as these dissipative effects are more pronounced in the nonlinear regime 
near the final coalescence, the separation between the various mechanisms is inevitably of an approximate nature. 
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1. tidal heating 



The first kind of dissipation we consider originates from the tidal deformation of m as the hole m falls in the static 
gravitational field generated by M. In the membrane paradigm |Q of black holes, in which the horizon is treated as 
a 2-D surface living in a 3-D space, endowed with physical properties like viscosity, this tidal deformation heats up 
the horizon. The heating is described by the horizon equations BJ-pq] 



dt 



CTab + [g - 0)aab + {2(7ac + la 



(25) 



(26) 



-Jilab ^'^Crab+ labk 

at 



(27) 



Here is the 2-D metric of the horizon of the infalling hole m, g = l/(4m) is the surface gravity of the hole m, and 
eab is the normalized electric part of the Weyl tensor {Cafj.bi^l^l'^ , with the horizon generators). 9 is the expansion 
rate of the horizon generators 



1 d 
'AAdt 



AA, 



and aab is the shear of the horizon generators 

d 



^^{AsY = 2{nab + -07afc)Ax"Ax^ 



(28) 



(29) 



with Ax'^ the coordinate separation of the horizon generators. 

For a hole with mass m falling in the external tidal field M/r^ , eab in an orthonormal basis can be approximated 

by [0 



M 



(30) 



where V is the velocity of the fiducial observers on the horizon of m moving in the external tidal field {V = ^j2M/r 
for infall from 00). Since <Jab and in Eqs. ( ^5| - |2^ ) are driven by e-g, which is small in our case, all nonlinear terms 
in Eqs. ( p5[p7| ) can be dropped, and "fab is decoupled. 9 is approximately given by |^ 



<Jab{t')a''\t')G{t,t') 

2M 



dt' 



<Jab<j'''G/r 



dr 



(31) 



(Jab = I eab{t')G{t, t')dt 
[eabG/r]dr , 



where G{t, t') is the teleological Green's function 

G(t,t') . 



exp[g(t - t')] for t < t' 
for i > t' 



(32) 



(33) 



is infinite where caustics exist on the horizon, although its contribution to the total increase of the horizon area is finite. 
We do not consider the effect of caustics in this paper. 
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In Eqs. ( |3l| ) and (|32|) , f is given by Eq. jH); the integrands are regarded as functions of r, with t = t{r) obtained by 
integrating Eq. (|24[). In Eqs. (|l]) and (|32|) the integrations are cut off when the hole falls through the horizon of M. 
In principle the integration should be carried over all times, but the residue is unimportant for the present purpose. 
In Fig. |l^, we plot a^^ (solid line) and 9 (dotted line) as functions of the Schwarzschild coordinate r, for the case of 
m = M. The horizontal and vertical scales are in terms of M = to = 1. 

A fraction of the initial gravitational potential energy fh — Ato/to is dissipated into the heating of the horizon of 
TO by this effect: 

Ato lAA 1 /■ , , 

2m+2M 

[e/r]dr (34) 

Notice that the integration in Eq. ( |34| ) is terminated at the point when the two holes are engulfed by a common 
horizon. It is irrelevant whether the object that has fallen in is a point mass or a black hole after that point. We 
approximated that point to be when the holes are separated by 2m + 2M in Schwarzschild coordinates, i.e., when 
the two holes are nearly touching. For m = M, the heating on the horizon of M is the same as that on to, hence the 
total fraction of energy going into heating of the horizons as the holes are falling in each others tidal field is given by 
2ffi. In Fig. |l^. the reduction factor for the energy available for wave generation, 

Fh^l' 2fh (35) 

is plotted against Tq, the initial separation of the holes, for the case of to = M. We see that Fh is decreasing with 
increasing initial separation, as a larger initial separation leads to a larger velocity and in turn larger e-^ in Eq. (|3^). 
As To oo, Ffi decreases to 0.86. For the range of covered in the numerical simulation, this effect reduces the 
gravitational wave output by about 10%. 




2. absorption of gravitational waves 

The second kind of dissipation arises from the fact that, unlike a point mass, a black hole has finite size. As it 
sweeps through the spacetime, it can reabsorb the gravitational wave already generated in the spacetime. 

The gravitational wave that a black hole can absorb depends on the frequency uj of the wave and the I — pole of 
the wave: 

wave absorbed = / ^ P*^^^ wave incident on m) x Ti{uj)du} . (36) 

Here T/ is the transmission coefficient for incident I pole waves as calculated in black hole perturbation theory, see 
e.g., 1^. We take the ^-pole wave incident on the infalling black hole m as the ^-pole wave in the spacetime times the 
cross-section of the hole to, as seen from the hole M, i.e., 

. , 7r(2m)2 1 ^2 

cross section ^- — -77 = 37 

~ 47r(2TO + 2Mf 4 A'P ^ ' 

The energy in the ^-pole wave in the spacetime is given e.g., in Ref. and is reproduced in Fig. ^ (dashed line) for 
Z = 2 as a function of lo. The vertical axis (for the dashed line) is dE/duj, all in scale of Af = 1. The transmission 
coefficient times 0.1 is also given as a function of lu (dotted line). Notice that the quasi- normal frequency is 
0.3737to^^. T/=2 is about 0.5 at this point, dropping to zero rapidly for a smaller uj. As the peak of dE/dio is at 
a smaller w, the product of T dE/du) is small, given by the solid line in Fig. ^ for the case of to = M. The area 
under the curve is found to be ~ 0.0012to. Comparing to the total gravitational wave energy in the I = 2 mode [ p8| , 
Ei^2 — 0.0092to, Eqs. ( |3^ ) and ( ^ lead to about 1% reabsorption of the I = 2 wave energy. 

Reabsorption of higher / modes can be estimated similarly. However, since the peak of dEi/du is always at an uj 
less than the corresponding quasi-normal frequency w;, whereas the transmission coefficients rise to larger than 0.5 
only for u) > ui, the reabsorption is always a small fraction of the corresponding component. As the I = 2 mode is 
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the dominant component making 90% of the total radiation, we see that this reabsorption effect can only lower the 
total energy output at a 1% level f\ 

Fabs « 99%. (38) 

We point out that there are other mechanisms causing the heating of the black hole horizons during the coalescence 
of black holes. Indeed, some of them are much larger than the effects considered above. For example, in the late 
stages of the coalescence, the infalling hole m is moving relativistically and is beaming gravitational waves in the 
forward direction. Such beaming causes large horizon heating p5| , p6[ | . The energy dissipated is an order of magnitude 
larger than the total gravitational wave output to infinity. However, such an effect is already implicitly included in 
the result of Hence no modification factor is needed to account for the process. The reduction in gravitational 
wave output due to the interference from the different parts of a body for a finite size object falling into a black hole 
has been estimated in Ref . Q . For the case of a two black hole collision, the effect is negligible. 



D. comparison with numerical results 

In Fig. |l^ we plot the final result of the analytic estimate for total gravitational wave energy output 

E = Fr^FhF^^,, X 0.0104^ (39) 

versus the initial separation between the two holes, for the case of m = M. It is represented by the dashed line. For 
comparison to the final result of Eq. we have also plotted two intermediate results as straight, solid lines. A simple 
application of the standard perturbative "DRPP calculation" (Eq. ^5]) overestimates the energy by roughly an order of 
magnitude. Replacing m by the reduced mass /i as in Eq. ( pO| ) gives the "Reduced Mass Calculation" , accounting for 
the finite mass of the "perturbing" black hole. The three correction factors i^r„, ^X, and Fabs together act to reduce 
the energy output further, leading to the curved, dashed line in Fig. ^ In view of the various approximations one 
has to make to obtain Eq. (|3^), the agreement between the analytic and the numerical results is remarkable, as the 
analytic results were obtained without prior knowledge of the numerical results, and vice versa. For L/M less than 
about 9 the analytic formula overestimates the actual energy output computed numerically. This is to be expected 
since for small enough separations the holes are initially engulfed by a common effective potential (in the sense of the 
potential in the black hole perturbation theory), or even a common event horizon. The approximation for colliding 
black holes in these cases is inappropriate. 

The connected circles show the maximum possible radiation output obtained by comparing the initial black hole 
masses estimated by the areas of the horizons (or a single horizon if the holes are close enough) to the total mass of 
the spacetime. For large separations this number approaches 29% as expected from the work of Hawking p9[ |. 

What physical understanding is gained from this semi-analytic exercise? The central message is that, as far as 
gravitational radiation is concerned, a black hole falling into another black hole is not much different from that of 
a point particle falling into a black hole. As the leading order approximation, the energy output can be described 
quite well by Eq. which is the test point particle result obtained in |^ modified by insight from the quadrupole 
formula. Wc found that the biggest effect due to black hole internal dynamics can be understood as the deformation 
of the horizons when the holes are falling in the tidal fields of one another. We see that the energy of about 10^'^ Af 
(for black holes of mass M) is dissipated by the viscosity of the horizon. Instead of radiating out to infinity, this 
portion of the initial potential energy is dissipated into the holes, increasing the horizon area. The smallness of this 
number again testifies that the horizon of a black hole is "stiff" , making a black hole rather like a point particle. At 
present we can only say the contribution from horizon heating is consistent with our numerical results to the level 
of accuracy. A future direct confirmation of this effect, especially for the case of non-equal mass black holes, will be 
particularly interesting. 

There are other effects that we considered, for example the reabsorption of radiation (also interference calculated 
in 1^) due to the finite size of the black hole. We found these effects to be negligible beyond the present level of 
accuracy, adding further weight to the understanding of a black hole behaving to a large extent as a point particle. 



^As the effect of reabsorption is small, various refinements of equations ( p6[ ) and (^^ are not meaningful, e.g., the cross- 
section is in fact a function of the separation between the two holes and the 1 = 2 mode of the two black hole system is 
different from the 1 = 2 mode of the infalling hole of mass m. 
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There are also other effects that we have thought about but do not know yet how to calculate analytically, for example 
the increase of the horizon area of the black hole due to caustics. However, the agreement between the numerical 
results and our present semi-analytic approximations suggests that these other effects most likely do not affect the 
total energy radiated significantly. 



V. CONCLUSIONS 



We have performed numerical and analytic calculations predicting the gravitational waveforms generated and total 
gravitational wave energy emitted when two equal mass black holes collide head-on. Waveforms for all cases studied 
show similar behavior: the normal modes of the final black hole are excited and account for most of the emitted 
signal. Both the ^ = 2 and £ — 4 waveforms are fit nicely by a superposition of the fundamental and first overtone of 
the black hole quasinormal modes. Although the fit to perturbative calculations of the £ = A waveform is quite good, 
the amplitude and precursor of this waveform are more sensitive to the various computational parameters than the 
more dominant £ = 2 waveform. 

The total energy radiated is on the order of 0.002M (where M is half the ADM mass of the spacetime), far below the 
estimate given by a simple application of the area theorem. The analytic study, appropriate for holes that are initially 
separate, confirms and elucidates the numerical results. We find the total energy radiated can be approximated quite 
well using the point particle result modified slightly to account for mass scaling, finite initial separation and internal 
dynamics of the black holes. Taken together, the analytic and numerical results indicate that even for holes that are 
initially infinitely separated, the total energy output will be the same order of magnitude. For throats that are close 
together. Price and PuUin ||38| have treated the evolution via gauge-invariant perturbation methods, regarding the 
system as a single, perturbed hole. They find remarkable agreement with our work for the £ = 2 waveforms and total 
energy radiated, independently confirming our results and providing a new method for evolving distorted black hole 
data sets. A detailed comparison of the numerical results presented here to semi-analytic results based on regimes 
where the throats are initially near or far from each other will be published elsewhere. 

The work presented in this paper is only a first step on the long path to computing the fully general evolution of 
three-dimensional, spiralling, coalescing black holes. We expect the more general case to be significantly more complex 
to compute. In axisymmetry our current calculations can be extended to include: boosted black hole collisions in 
which the holes are given initial finite velocities, unequal mass black holes which can radiate not only gravitational 
waves but net linear momentum, and spinning colliding black holes where one can expect more energy to be radiated, 
particularly if the holes have opposite spin vectors. We intend to pursue such extensions to our present axisymmetric 
code and to develop and apply new, more general, three dimensional codes to these systems as well. 

Finally, we note that we have prepared a video showing results for several of the simulations reported here. Interested 
readers may contact NCSA media services at the internet address media@ncsa.uiuc.edu for information on how to 
obtain a copy of the video entitled "The Collision of Two Black Holes." At this address one can also obtain a copy of 
a videotape of the original movies based on the work of Smarr and Epp ley. The NCSA group has also set up a World 
Wide Web server accessible at the URL littp : / /j ean-luc . ncsa. uiuc . edu , and there one can find images and movies of 
the simulations presented in this paper that cannot be published in traditional form. 
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FIG. 1. A single quadrant of the Cadez grid is displayed for the case /x = 2.2. The throats are centered on the axis at 
z = ± coth fi. Lines of constant rj concentrically surround the throats locally, and become spherical far from the holes. 

FIG. 2. The metric function 7,,,, — A for the case /i = 1.2 is shown at coordinate time t = 25M, where M is the mass 
parameter defined to be half the ADM mass. This configuration was surrounded by a global apparent horizon initially. 

FIG. 3. The lapse function a is shown for the case /x = 1.2 at the time t = 25M. It collapses to zero uniformly around the 
entire region surrounding both throats, indicating that the system is behaving as a single black hole. 

FIG. 4. The lapse function a is shown for the case fj, = 3.25 at time t = 22. 5M. At this point the two black holes are still 
evolving independently of each other. 

FIG. 5. The mass of the apparent horizon is shown for the cases jj, = 1.2 (solid line) and jj, — 2.2 (dashed line). In the 
H = 1.2 case the initial horizon surrounds both holes and most of the mass-energy is contained within the horizon. After a 
slight increase early in the evolution (within the first few M) the mass of the horizon remains essentially constant until t ~ 30M 
when it begins to grow due to numerical effects discussed in the text. In the fj, = 2.2 case the initial horizons around each 
throat are distinct. At t ~ 17M a new apparent horizon appears, surrounding both holes and accounting for essentially all the 
mass-energy in the system. 

FIG. 6. The area of the constant 77 surface that just encircles the saddle point is plotted as a function of ix. The solid line is 
IGtt, the area corresponding to the total ADM mass of the spacetime (we work with units such that the total mass is normalized 
to unity). 

FIG. 7. We plot various timings of the evolution for the six initial data sets, tnewt is the Newtonian free fall time required 

for two particles to collide from rest at this separation, thor is the time at which the apparent horizons merged, and tring is the 
time at which the first (negative) peak is seen in the i = 2 Zerilli function recorded at r = 40M. 

FIG. 8. The £ = 2 waveforms for the case n = 1.2. The solid line is the numerically generated waveform extracted at 
r — 40M. The dashed line is a fit of the two lowest i = 2 quasinormal modes, over the domain 70 < t/M < 160, to the 
extracted waveform. 

FIG. 9. The £ = 2 waveforms for the case fj, = 2.2. The solid line is the waveform extracted at r = 40M and the long 
dashed line is the waveform at r = 60M. The dotted line is the quasinormal mode fit. 

FIG. 10. The £ — 4 waveforms for the case /i = 2.2. The solid line is the waveform extracted at r = 40M and the short 
dashed line is the quasinormal mode fit. 

FIG. 11. The £ = 2 waveforms for the case fj, = 2.7. The solid line is the waveform extracted at r = 40M and the short 
dashed line is the quasinormal mode fit. 

FIG. 12. The £ = 4 waveforms for the case n = 2.7. The solid line is the waveform extracted at r = 40M and the short 
dashed line is the quasinormal mode fit. 

FIG. 13. The time evolutions for three different radiation indicators are shown for the case // = 2.2. The solid line is the 
Zerilli function (£ = 2), the dashed line is the Newman-Penrose scalar ^'4, and the dotted line is the radial component of the 
Bel-Robinson vector. The data tracks the behavior in time of a single point on the grid located at r = 70M along the equator. 
The £ = 2 fundamental quasinormal mode is clearly present and dominant in all three signals. 
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FIG. 14. The total gravitational wave energy output is shown for the six parameter studies. The connected circles are the 
upper limit based on the area theorem, the clustered symbols show numerical results at various detector locations, and the 
crosses show early results by Smarr and Epplcy with their approximate error bars. The sohd Hnc labelled "DRPP Calculation" 
is the result of a naive application of the point particle result, the solid line labelled "Reduced Mass Correction", takes into 
account the finite mass of the infalling black hole, and the dashed line is the semi-analytic estimate, including several effects 
discussed in the text. Finally, the dot-dashed line is an empirical fit to the data for low values of the separation parameter /i, 
showing the exponential falloff of the energy in the one black hole regime. 

FIG. 15. We plot the energy factors Fr^ and Fh vs. Vo- The horizontal axis is the initial Schwarzschild coordinate To of the 
infalling hole in units of M. Fr„ is the gravitational wave energy output for infalling from ro divided by that from oo. Fr„ is 
smaller for small To mainly because the velocity of the infall is smaller. Fh is the percentage of the available energy compared 
to the total gravitational energy available after taking into account the dissipation due to tidal distortions of the black hole 
horizons. As To — > oo, F/, tends to 86%. The curves are calculated for m = M. 

FIG. 16. Development of horizon shear and expansion. The horizontal ax;is is the Schwarzschild coordinate r of the infalling 

object in units of M . The solid line represents the magnitude of the horizon shear on the horizon of the infalling hole m. The 
dotted line is the expansion. Both the shear and the expansion are in units of and are plotted for the case m = M. 

FIG. 17. The reabsorption of the I — 2 wave. The horizontal axis is the angular frequency cu of the / = 2 wave in units 
of M^^ . The dashed line is dE/dio, the dotted line is the transmission coefficient of the I = 2 wave multiplied by 0.1 and the 
solid line is the product T{dE/dio). For M = m the area under the solid curve, which is found to be 0.0012m, is roughly the 
total reabsorption of the I = 2 wave by the infalling hole. 



M 


M 


L/M 


Apparent horizon 


Event horizon 


1.2 


1.85 


4A6 


global 


global 


1.8 


0.81 


6.76 


separate 


critical 


2.2 


0.50 


8.92 


separate 


separate 


2.7 


0.29 


12.7 


separate 


separate 


3.0 


0.21 


15.8 


separate 


separate 


3.25 


0.16 


19.1 


separate 


separate 



TABLE I. The physical parameters for the six initial data sets. M is the mass parameter equal to half the ADM mass of 
the spacctime, L/M is the proper distance between the throats, and we note whether or not a single apparent or event horizon 
surrounds both holes. 





E^ 


£*4 


Epy 


1.2 


1.34 - 0.99 xlO"" 


9.81 - 1.32 xlO"*' 


7.58 - 1.29 xlO"'' 


1.8 


1.95 - 1.64 xlO'" 


3.46 - 1.72 xlO"" 


3.77 - 1.73 xlO"* 


2.2 


6.10 - 5.27 xlO"" 


7.76 - 5.21 xlO"" 


8.91 - 5.26 xlO"* 


2.7 


7.50 - 6.86 xlO"" 


7.16 - 5.00 xlO"* 


7.52 - 5.07 xlO"* 


3.0 


8.85 - 7.13 xlO"" 


9.22 - 4.86 xlO"" 


10.0 - 4.93 xlO"* 


3.25 


1.37 - 0.85 xlO"-' 


4.32 - 1.11 xlO"'' 


3.39 - 1.11 xlO"^ 



TABLE II. The total radiated energy for the six initial data sets normalized to the ADM mass (2Af ) of the spacetime. We 
compare results using several different methods of calculation on a grid resolved with 200x35 zones. The energies are extracted 
at five different radii, and the range given here is from the innermost to the outermost radius. See text for details. 
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